Searching for vector bileptons 
in polarized Bhabha scattering 
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In this paper we analyze the effects of virtual vector bileptons in polarized Bhabha scattering at 
the energies of the future linear colliders. In order to make the calculations of the differential cross 
sections more realistic, important beam effects such as initial state radiation, beamstrahlung, beam 
energy and polarization spreads are accounted for. The finite resolution of a typical electromagnetic 
calorimeter planned for the new linear colliders is also considered in the simulation. The 95% 
confidence level limits for bilepton masses in 331 models are evaluated. 
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I. INTRODUCTION: 

Although the standard model (SM) explains all current 
experimental data, it is commonly believed that it is not 
the final answer. Many questions have been left open, 
so many theorists believe that there should be some new 
physics lurking at the TeV scale. Many models have been 
proposed to explain these unanswered questions. Ex- 
amples are supersymmetry, composite models and grand 
unification theories. All extended models have in com- 
mon the prediction of new particles. This paper concerns 
nonstandard particles known as bileptons |l[ - bosons 
which couple to two leptons and which carry two units of 
lcpton number. They are present in several models, such 
as technicolor, left-right symmetric models, and in the 
grand unification schemes. Heavy gauge bileptons appear 
in extended gauge models where the electroweak group 
is imbedded in a larger group. The so-called 331 models, 
based on the gauge group SU(3)c O SU(2)l ® C/(l)y, 
fall in this category. In this paper we discuss how to 
search for signals of an off-mass-shell, doubly-charged 
vector bilepton in polarized Bhabha scattering, at the 
new linear collider energies. Most of the analysis is made 
in the framework of 331 gauge models Q, but the pre- 
dictions of an SU(15) model [3(] are also presented, and 
may be taken as a measure of the mo del- dependence of 
our results. The detection of such indirect effects of a 
vector bilepton would be strong evidence of new physics. 
Should any signal of a vector bilepton be detected in a 
high energy collider, it would be important to settle the 
question of the correct underlying model. 
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For realistic comparison with experiment, we consider 
important beam effects such as (i) initial-state radiation, 
beamstrahlung and beam energy spread; (ii) longitudi- 
nal and transverse polarization of the colliding beams, 
which are expected to be available in the next gener- 
ation of linear colliders; (iii) Gaussian smearing of the 
four-momenta of the final-state leptons, simulating the 
uncertainties in the energy measurements in the electro- 
magnetic calorimeters. In section II we give a brief re- 
view of the 331 minimal model |4j. Section III describes 
in detail the numerical simulation of the Bhabha events. 
In section IV we analyze several observables such as an- 
gular distributions and asymmetries, with the purpose 
to look for indirect vector bileptons signals. In section V 
we establish bounds on the couplings and masses of these 
bileptons at the 95% confidence level, from the angular 
distribution of the final-state leptons. As we are working 
in the context of the minimal 331 version, one can also 
verify the validity of the relations between the masses of 
the vector bileptons and new neutral gauge bosons, which 
are connected to the Higgs structure of the 331 models. 



II. REVIEW OF THE 331 MINIMAL MODEL 

The fcrmionic content of the minimal 331 model ar- 
ranges the ordinary leptons in SU(3)l antitriplets, two 
generations of quarks in triplets and the third genera- 
tion in an antitriplet. The anomaly cancellation takes 
places only when all three families of quarks and leptons 
are taken together. The model also predicts a new neu- 
tral gauge boson Z' and four vector bileptons Y ± and 
V ±=1= , which acquire masses after spontaneous symme- 
try breaking (SSB): first SU(3) L x U(1) Y breaks down 
to SU(2)l x U(1)y- This can be accomplished with 
only one SU(3)l scalar triplet; the next stage of SSB, 
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SU(2)l x U(1)y — > U(l)c, requires an additional scalar 
SU(3)l triplet. In order to provide quarks and leptons 
with acceptable masses, however, two additional triplets 
and a sextet would be necessary for symmetry breaking 
with a minimal Higgs structure. In this minimal version 
of the 331 model, the mass My of the doubly-charged 
bilepton Y ++ and the mass Mz* of the neutral vector 
boson Z' are correlated parameters given by the expres- 
sion: 

M Y ^3(1 - 4 sin 2 6 W ) 
M z - ~ 2cos6>w 

This relation no longer holds for the SU(lb) or any 
other model with different Higgs structure. Presently, 
the most useful lower bounds on the vector bilepton 
mass can be derived from fermion pair production at 
LEP and lepton-flavor violating charged lepton decays 
5], M Y > 740GeV, and from muonium-antimuonium 
conversion [||, My > 850GeU. While the former is less 
stringent, it does not depend on the assumption that the 
bilepton coupling is flavor diagonal. The data collected 
by the CDF Collaboration at the Fermilab Tevatron ex- 
clude 331 Z' masses below 920 GeV 0]. 

It has been argued [8( that the minimal 331 model 
described above can not be analyzed with perturbation 
theory around an energy scale /i, at which the ratio 
9x1 'A = sin 2 9 w (fi)/(l - Asin 2 6 w {p)), where g x and 
Qh denote the U(l)x and SU(3)l gauge group coupling 
constants, develops a Landau-like pole. The scale [i was 
found to be of the order of 4TeV for the minimal 331 
model. Solutions |9J have been put forward to circunvent 
this possible nonperturbative nature of some 331 models 
at the TeV energy scale. In particular, by adding three 
octets of vector leptons to the particle content, it is possi- 
ble to make sin 2 9w{n) decrease with the increase of the 
energy scale as far as the light states of the standard 
model are concerned, and thus the perturbative regime 
is restored for the TeV energy scale. 

III. SIMULATION 

When a charged particle is accelerated to energies 
much superior than its own rest energy, the probability 
of photon emission is very high. The so called initial- 
state radiation (ISR) is the emission of photons by the 
incoming electrons and positrons due this acceleration. 
ISR is a dominant QED correction in leptonic scatter- 
ing. Using the structure function approach discussed in 
10], ISR effects were taken into account in this simula- 
tion. We neglect electroweak loop corrections throughout 
the calculations. For the energies planned for the future 
linear colliders, the first-order weak loop-corrections to 
the Bhabha angular distribution [ll[ are at most of the 
order of 10% of the corresponding distribution at Born 
level, for wide-angle scattering. The remaining weak 



loop-corrections are smaller. We expect that the mag- 
nitude of weak loop-corrections in the extended models 
we are studying be similar to that of the standard model. 
Moreover, to some extent these corrections should cancel 
out in the asymmetry ratios, and therefore their effect on 
the observables discussed in the next section should be 
small. In the quest for higher luminosities at the linear 
colliders, very dense incoming beams, with transverse di- 
mensions of the order of a few dozen nm, will be needed. 
The result is a large charge density, leading to very strong 
electromagnetic fields. For this reason, the particles in 
a colliding bunch suffer considerable transverse accelera- 
tion, which gives rise to the emission of synchrotron radi- 
ation, the so-called beamstrahlung. The effective energy 
available for the reaction is then smaller than the nominal 
value. Machines with large luminosity per bunch crossing 
produce more beamstrahlung. Thus, the average energy 
loss of a positron or electron by beamstrahlung depends 
on the design parameters of the accelerator. For some 
designs the beamstrahlung energy loss may even reach 
30% of the available energy [l2j • Using the approach dis- 
cussed on ref. [l3| we obtained the beamstrahlung struc- 
ture function, starting from an energy-dependent set of 
NLC design parameters [l4 |. Convoluting the ISR and 
beamstrahlung emissions spectra we obtained the distri- 
bution which was used to compute the required differ- 
ential cross sections. Our simulations also included a 
Gaussian-distributed beam energy spread, with a width 
of 1% of the nominal beam energy. 

As we are working mostly in the framework of the 331 
minimal model, two new Feynman diagrams have to be 
added to the standard diagrams for Bhabha scattering: 
a s-channel exchange of a Z' boson and a u-channel ex- 
change of a doubly-charged bilepton. As a partial check 
of our calculation, we "switched off" the beam polar- 
izations, ISR and beamstrahlung effects and the Z' ex- 
changes as well, and verified that we indeed reproduce 
the results of the trace calculation of ref. [l5j]. Like- 
wise, we cross-checked our calculations with those of ref. 
[l6| , which were carried out in the framework of the stan- 
dard model, but with arbitrary beam polarization. Us- 
ing Monte Carlo techniques we calculated the differential 
cross sections with the simulated events selected accord- 
ing to the following set of cuts : (i) the final-state elec- 
trons and positrons were required to be produced within 
the angular range |cos6*i| < 0.95, where 9i stands for the 
polar angle of the final-state lepton with respect to the 
direction of the incoming electron beam; (ii) all events in 
which the acollinearity angle £ of the final-state e + — e~ 
three-momenta did not pass the cut C < 10° were re- 
jected; (iii) the ratio of the effective center-of-mass en- 
ergy to the nominal center-of-mass energy for any ac- 
ceptable event was required to be greater than 0.9. We 
considered a 500 fb~ x integrated luminosity, which could 
be achieved in several years of machine operation , con- 
sidering a maximum center-of-mass energy of 1 TeV. 
Finally, we simulated the finite resolution of the NLC 
electromagnetic calorimeters by Gaussian-smearing the 
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four-momenta of the produced electrons and positrons 
[l7l ]. The directions of the lepton three- momenta were 
smeared in a cone around the corresponding original di- 
rections, whose half-angle is a Gaussian with half-width 
equal to 10 mrad, whereas the energies of the final-state 
lcptons were distributed as a Gaussian with half-width 
AE of the form AE/E = 10%/VE © 1%. 

As the experiments at SLC have shown, beam polar- 
ization is an important tool in the search for new physics 
and for precision measurements in particle physics. Elec- 
tron and positron beams with high degrees of polariza- 
tion will play a significant role at the new linear collid- 
ers [l8j]. Although lepton polarization has the inconve- 
nience of being accompanied by considerable beam in- 
tensity loss, the new linear colliders will compensate by 
achieving very high luminosities. Polarized lepton beams 
can reduce backgrounds and increase of the sensitivity of 
spin-dependent observables to nonstandard phenomena. 
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FIG. 1: Total cross section as a function of the center-of- 
mass energy ^/s, for My = 1.2TeV (331); solid line (SM) 
represents the standard model cross section. 



IV. OBSERVABLES 



In our numerical calculations, the longitudinal 
(transversal) beam polarizations were taken to be Pl = 
90% (P T = 90%) for electrons and P L = 60% 
(Pt = 60%) for positrons, with an uncertainty given by 
AP L /P L = 0.5% (APt/Pt = 0.5%). Firstly, we only 
considered longitudinal polarization of the beams. 

FigQ]shows the dependence of the total cross section of 
the 331 minimal model on the center-of-mass energy, for 
an input bilepton mass My = 1.2TeV. For acceptable 
values of the bilepton mass, the deviation from the stan- 
dard model cross section is small over the energy range 
expected for the next generation of linear colliders. This 
was also found to be true for the SU (15) model. 

Next we looked at the angular distributions in the min- 
imal 331 model at y/s = 500GeV and ^/s = ITeV, which 
are displayed in Fig[^]with the corresponding curves for 
the standard model and SU (15) shown for comparison. 
All curves indicate that the final-state electrons are emit- 
ted mostly in small angles (with respect to the incident 
electron beam), but significant differences among the an- 
gular distributions still occur over the remaining angular 
range. These deviations from the SM predictions could 
be helpful to establish the correct underlying electroweak 
model, should any signal of a vector bilepton be detected 
in a high energy collider. In a previous work [2C| . a sim- 
ilar pattern was found for the sensitivity of M0ller scat- 
tering total cross section and angular distribution to the 
coupling of vector bileptons. 

As the angular distribution for 331 and standard mod- 
els are strongly asymmetric, a typical non-identical parti- 
cle final state behavior, the integrated forward-backward 
asymmetry is large: 
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FIG. 2: Angular distribution of the final-state electrons; 
Curves (SM-5) and (SM-10) show the angular spectra pre- 
dicted by the standard model at = 500GeV and \fs = 
ITeV respectively, while (331-5) and (331-10) display the cor- 
responding angular distributions in the minimal 331 model. 



where z stands for the cosine of the angle 9 between the 
outgoing and incoming electrons. This is displayed in 
FigGH where Afb is plotted for several bilepton masses 
at an energy \fs — 500Gey. As expected the 331 Afb 
value approaches the standard model Afb value as the 
bilepton mass increases. 

Also, we analyze the discovery potential of spin asym- 
metries. 



J° dz— — J 95 dz— Using the polarization-dependent angular distribu- 

Apb — n°' 95 — tttt^ — (1) tions, we define the following asymmetries: 

J-0.95 ""^ dz ~ JO a ^ dz 



4 



A^cosO) 



da(- 



\P b L \) + da(-\Pl\, \P b L \) - da(\P£\, -\P b \) - da(\P£\ 



da{- 



\P£l- 



\P b L \ 



da{- 



M\Pl 



1^1 



1^1 



■M\p£\>\n\) 



(2) 



1000 
M Y (GeV) 




FIG. 3: Forward-backward asymmetry Afb for several input 
masses My, at ^/s — ITeV, according to the minimal 331 
model (331). Curve (SM) shows the standard model predic- 
tion 
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FIG. 4: Polar angle dependence of spin asymmetry Ai(cos8) 
at y^s = ITeV; Curve (SM-10) shows the standard model 
predictions for A\{cos6), while (331-10) represents the corre- 
sponding expectations for the minimal 331 model, for a mass 
M Y = 1.2TeV. Curve (SU15-10) shows the prediction for the 
SU(15) model. 
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FIG. 5: Polar angle dependence of spin asymmetry A2(cos9) 
at y/s — 500GeV; Curve (SM-5) shows the standard model 
predictions for A2(cos8), while (331-5) represents the corre- 
sponding expectations for the minimal 331 model, for a mass 
My = 1.2TeV. Curve (SU15-5) shows the prediction for the 
SU(15) model. 



In the definitions of the three asymmetries above, a 
and b refer to the electron and positron beams respec- 
tively. Fig 0] shows the angular dependence of asymmetry 
Ai for a bileptonic mass My = 1.2 TeV, at -y/s = ITeV. 
We do not show the corresponding 500 GeV curve, since 
the Ai sensitivity is considerably lower at this energy. 
Asymmetry Ai has proved to be quite sensitive to the 
presence of vector bileptons. For y/s = 500GeV this is 
shown in Fig|5j where the A2 dependence on cos 9 is seen 
to be weaker in the SM, than in the two other models 
endowed with bileptons. The discriminating power of 
Ai is even more evident at y/s — ITeV , as depicted in 
FigEl The resulting errors associated with the calcula- 
tion of longitudinal asymmetries Ai and A 2 were small, 
and hence the error bars are not shown in the figures. In 
particular, using bootstrap resampling techniques (l9| we 
confirmed that statistical errors are indeed small. Some 
of the systematic errors, such as the uncertainties on the 
degrees of beam polarization, were directly included in 
the simulation. As a matter of fact, we expect some of 
these systematic effects to cancel out in the polarization 
asymmetries. 

Asymmetry A3 concerns the difference between a po- 
larized distribution and its unpolarized counterpart, and 
is displayed in Figsl7]and[8]for both 1 TeV and 500 GeV 
respectively. For both energies its the 331 minimal model 
to give the most different predictions compared with the 
SM. A curious fact that must be noted it's the change of 
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FIG. 6: Polar angle dependence of spin asymmetry A2(cos9) 
at — ITeV; Curve (SM-10) shows the standard model 
predictions for A2(cos9), while (331-10) represents the corre- 
sponding expectations for the minimal 331 model, for a mass 
M Y = 1.2TeV. Curve (SU15-10) shows the prediction for the 
SU(15) model. 
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FIG. 8: Polar angle dependence of spin asymmetry A3 (cos8); 
Curve (SM-5) shows the standard model predictions for 
As(cos8) at t/s = 500GeV , while (331-5) represents the corre- 
sponding expectations for the minimal 331 model, for a mass 
M Y = l.2TeV. Curve (SU15-5) shows the prediction for the 
SU(15) model. 
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FIG. 7: Polar angle dependence of spin asymmetry Az(cos9); 
Curve (SM-10) shows the standard model predictions for 
As(cosO) at y/s = ITeV, while (331-10) represents the corre- 
sponding expectations for the minimal 331 model, for a mass 
M Y = 1.2TeV. Curve (SU15-10) shows the prediction for the 
SU(15) model. 
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FIG. 9: Transverse polarization asymmetry Ait(cos8) as a 
function of cos8, at ^/s = ITeV. The line with error bars 
(SM) corresponds to the standard model prediction, while 
the other curves represent the asymmetry for the minimal 
331 model and SU(15) model with M Y = 800GeV. The error 
bars for (SU15) and (331) are similar to those of the SM on 
a bin-to-bin basis. 



behavior that only the 331 model suffers at 1 TeV . At 
500 GeV the three curves have qualitative similar distri- 
butions, but at 1 TeV, although this similarity holds for 
the SM and the SU(15) model, this is no longer valid for 
the 331 model, as one can see in FigJTJ 

Finally, we analyze the discovery potential of the spin 
asymmetries Ait and A^t for Bhabha scattering. We 
define: 



. , n s ^ d(cos8)d<t> /(-) °^ d(cos0)d<p 
AiT(COSO) = — tj- — ' 

J(+) ^ d(cos8)d<t> + J(-) °^ d{cos8)d<t> 

(5) 

where the subscript +(— ) indicates that the integration 
over the azimuthal angle <j) is to be carried out over the 
region of phase space where cos 2(f> is positive (negative) . 

Fig[5] shows the result for Ait considering My = 
800 GeV and energy of 1 TeV. Although the errors are 
bigger than the case with longitudinal polarization, the 
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FIG. 10: Transverse polarization asymmetry Ait as a func- 
tion of My at ITeV. The solid histogram (331) refers to 
the minimal 331 model. The resulting curve for the SU(15) 
model (SU15) is also shown for comparison. 
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FIG. 11: 95% C. L. contour plots on the (j 3 i/e,tfy) plane 
for longitudinally polarized Bhabha scattering, at the NLC 
center-of-mass energies y/s — 500GeV (lower curve) and y/s = 
ITeV (upper curve). 



observable allows a clear vector bilepton signal for the 
considered extended models. The changes in the shapes 
of the curves for A\t(cos0) for the two nonstandard 
models are small when the bilepton mass is increased 
to 1.2 TeV. 

The integrated version of the asymmetry A2T, were 
also investigated 

_ J {+) d(eos6)dcp - / ( _ } dicoaOW 

-r^2T — ~j2 7l2 

/(+) d{cos9)d^ + d(coa9)M 

(6) 

where the integrations are consistent with the cuts spec- 
ified in section III. The result is shown in FigJTOj where 
the mass dependence of An for \fs = ITeV for both the 
£?7(15) and the 331 model were plotted. Qualitatively, 
there are no differences between the curves. Nevertheless 
the absolute value of A^t is greater for the SU (15) model. 
For the sake of comparison, the approximate value of A^t 
for the standard model was found to be —0.1 . 



V. BOUNDS 

We estimated lower bounds on the bilepton mass by a 
X 2 test, comparing the angular distribution do / d(cosO) 
of the final-state electrons, modified by the presence of a 
vector bilepton, with the corresponding standard model 
distribution. Here only longitudinal beam polarization 
was considered. We assumed that the experimental data 
will be well described by the standard model predictions. 
We treat 173; and My as free parameters in an effective 
theory and define our x 2 estimator as 

^5 / N SM at \ 2 

X^,M Y )=^\^) (7) 



where is the number of standard model events de- 

tected in the i th bin, Ni is the number of events in the 
i th bin as predicted by the model with bileptons, and 

ANf M = \J(\/N? M ) 2 + (N? M e) 2 the corresponding 
total error, which combines in quadrature the Poisson- 
distributed statistical error with the systematic error, 
which for the latter we assumed to be of e = 5% for 
each measurement. The angular range |cos#| < 0.95 was 
divided into = 20 equal- width bins. For fixed values 
o'f the coupling g^i the bilepton mass My was varied as 
a free parameter to determine the \ 2 distribution. The 
95% confidence level bound corresponds to an increase 
of the x 2 by 3.84 with respect to the minimum x 2 m in °f 
the distribution. FigQT] presents the resulting 95% C. L. 
contour plots on the (gzil 'e,My) plane for the nominal 
center-of-mass energies y/s = 500GeT^ and y/s = ITeV. 
The unpolarized case is shown in Figd^l We also calcu- 
lated the corresponding 95% C. L. limits on the bilepton 
mass at these ILC energies, considering only the minimal 
331 model. The results are displayed in Table I. 

VI. CONCLUSIONS 

We found that the analysis of spin asymmetries in 
Bhabha scattering leads to tighter bounds than the coun- 
terparts obtained from Miller scattering [20| . We em- 
phasize that several beam effects such as ISR, beam- 
strahlung, beam energy and polarization spreads were 
taken into effect in our simulations. In brief, we showed 
that polarization-dependent observables may be used to 
search for signals of vector bileptons in Bhabha scatter- 
ing, at the next generation of linear colliders. Should 
vector bileptons be detected, the underlying nonstandard 
model describing their interactions could be determined 
by the discriminating power of combined measurements 
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TABLE I: 95% C. L. limits on the bilepton mass in the min- 
imal 331 model, at NLC energies 
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FIG. 12: 95% C. L. contour plots on the {gu/e, My) plane 
for unpolarized Bhabha scattering, at the NLC center-of-mass 
energies \fs = 500GeV (lower curve) and yfs = ITeV (upper 
curve) . 



of spin asymmetries. Otherwise, polarized Bhabha scat- 
tering would still allow that useful experimental bounds 
be set on the bilepton mass and couplings. 
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